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We derive the planar large N non-supersynimetric background of 
K> I the quantum mechanical hamiltonian of two hermitean matrices cou- 

j^ ■ pled via a Yang-Mills interaction, in terms of the density of eigenvalues 

of one of the matrices. This background satisfies an implicit non linear 
integral equation, with a perturbative small coupling expansion and 
a solvable large coupling solution, which is obtained. The energy of 
system and the expectation value of several correlators are obtained 
in this strong coupling limit. They are free of infrared divergences. 
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1 Introduction 

The study of multi-matrix model^j, and particularly their large N limit [1] , 
is of great interest. It is well known, for instance, that the large N limit of 
their description of DO branes [2] has been conjectured to provide a definition 
of M theory |3j. In the context of the AdS/CFT duality [1], [5], [6], due to 
supersymmetry and conformal invariance, correlators of supergravity and 
1/2 BPS states reduce to calculation of free matrix model overlaps [7j, ^ or 
consideration of related matrix hamiltonians [9]. For stringy states, in the 
context of the BMN limit [10] and A/" = 4 SYM, similar considerations apply 
HH, in], US]- A plane- wave matrix theory [H] is related to the A/" = 4 SYM 
dilatation operator [TJ]. Recently, multi-matrix, multi-trace operators with 
diagonal free two point functions have been identified [16], ^Y]. In earlier 
works [18], [19], [20] it has been argued that QCD can be reduced to a finite 
number of matrices with quenched momenta. 

In this communication we will consider the quantum mechanics of two 
hermitean matrices with harmonic potentials, interacting through a standard 
Yang-Mills quartic potential. 

We will use the approach, first developed in [21], of treating one of the 
matrices, which generates the large N planar background, in a coordinate 
representation, and the other in a creation/anihilation basis. This was done 
in the context of 1/2 BPS states and a dual free harmonic Hamiltonian [S], 
[9] with the matrix generating the background being the holomorphic com- 
ponent of a complex matrix. A precise phase space identification between 
the collective density description of the dynamics of this matrix [22], and the 
droplet description of the LLM [23] metric is obtained [2T]. The generaliza- 
tion of this approach to include qym interactions was developed in [23]. By 
considering the planar background generated by of one of the two Hermitean 
matrices, further properties of the spectrum were established in [25] . 

In [21], [23] and [25], a supersymmtric setting was always assumed, al- 
lowing one to consistently neglect normal ordering terms. As a result, the 
planar background is harmonic, and qym independent. 

In this communication we explore the consequences of not requiring super- 
symmetry, and establish a gvM dependent density description of the planar 
large N limit of the system. The communication is organized as follows: in 



^By matrix models we m.ean integrals over matrices or the quantum mechanics of matrix 
valued degrees of freedom 



Section 2, the hamiltonian of the two interacting harmonic matrices is intro- 
duced. It is shown how a Bogohubov transformation brings the hamiltonian 
to a form quadratic in the creation/annhihation oscillators of one of the ma- 
trices, with frequencies having a square root dependence on the eigenvalues 
of the other matrix [24j q The sector of the Hamiltonian contributing to the 
planar large A^ background is identified. In Section 3, the planar background, 
in terms of the density of one the matrices, is obtained implicitly through 
a non-linear integral equation. In Section 4, the perturbative expansion in 
A = Qym is described, and the first order correction obtained and shown to 
agree with perturbation theory. In section 5, the strong coupling limit of the 
background is obtained explicitly. In Section 6, it is argued that this strong 
coupling limit is indeed the background for two free hermitean matrices in- 
teracting through the Yang-Mills potential. This solution is free of infrared 
divergences. Section 7 is reserved for a summary and discussions. 

2 Model Hamiltonian and planar large N sec- 
tor 

We consider the quantum mechanical hamiltonian 



1 ?// 1 in 

H ^ -Tr(P,2) + -Tr(X,2)+-Tr(P|)+-Tr(X,2)-^^^,Tr([Xi,X2][Xi,X2]), 

(1) 
where Xi and X2 are two N x N hermitean matrices, and Pi and P2 their 
conjugate momenta, respectively. 

One can think of ([1]) as associated with two of the six Higgs scalars of the 
bosonic sector of A/" = 4 SYM, in the leading Kaluza Klein compactification 
on 5*^ X R. The harmonic potential results from the coupling to the curvature 
of the manifold. It should be borne in mind that we do not require super- 
symmetry. Alternatively, compactification of QCD2+1 on a sphere results in 
a similar hamiltonian. 

We will follow the approach first suggested in pT] of treating one of 
matrices, Xi, in coordinate space and exactly (in the large N limit), and the 
other, X2, in a creation/annihilation basis. Letting 



^See also |26l 



X,^^{A, + A\) P, = -^ l^{A, - Al) (2) 

the hainiltonian ([1]) takes the form 



_ 3-hLTr{2[X,,Al][X,,A2] + [X,,A,]' + [X,,Al]') (3) 

w w 

As the interaction is quadratic in the oscillators, one can perform a Bogoli- 
ubov transformation 

{V^A2V),, = cosh(0,,)5y- - sinh(0,,)4 (4) 

with 

tanh(20i,) = ^4^^ ^- , (5) 

where the Aj's are the eigenvalues of the matrix Xi and V is the unitary 
matrix that diagonalizes Xi. Then (|3]) takes the form 



1 7/;2 ^ / 1 

^=2Tr(A') + YTK^i)+EV^' + 2^^M(A.-A,)2 (4^..+ 2)- (6) 

We are interested in the leading large A^ configuration of the system. In the 
B,B^ sector of the theory, states with non zero B quanta lead to excited 
states and spectra, which contributions which are subleading in -^. The 
only contribution to the large A^ ground state configuration comes from the 
zero point energies of the B,B^ oscillators, and we are therefore led to the 
Hamiltonian: 



1 2 1 ^ 

Ho = ^MP^) + Y^r{Xf) + -Y, ^w^ + 2g'y,,{X,-\,y (7) 



3 A self consistent non linear integral equa- 
tion for the large N background 

The Hamiltonian (jTj) is describes the dynamics of a single hermitean matrix, 
and the large N background can be described in terms of the density of 
eigenvalues, 

i 

as the minimum of the cubic collective field effective potential [22] 

Veff ~ ~^ dx(l)^{x) + -^ dx(l){x)x^ ~ f^{ dx(l){x) — N) 

+ ^ dx dy^w^ + 2g^^j{x - yf 0(x)0(|/), (8) 

where the Lagrange multiplier /x enforces the contraint j dx(f){x) = N. To 
exhibit explicitly the A^ dependence, we rescale 



X -^ VN (f){x) -^ VN(j){x) 1^-^ N/j (9) 

and obtain 



Veff = N' 



7]-2 r ^2 

^ — / dx(j)'^{x) H / dx(f){x)x'^ ~ /^( / dx(j){x) — 1) 



+ \jdxj dy^w^ + 2\{x - y)20(x)0(y)] (10) 

where A = Oym^ is the usual 't Hooft's coupling. 

As A^ — > oo, the large N background configuration minimizes (ITOl) and it 
satisfies: 



vrVoW = 2^ - w'^x^ - 2 / dy^w'^ + 2A(x - yYMv) (H) 



4 Perturbative expansion 

When A = 0, (TTT!) reduces to the well known Wigner distribution: 



7l4>o{x) 



V 2u — 2u7 — w^x^ = V2w — w'^x^, \x\ < xq = \ — 

V w 



:i2) 



with the identification /i = 2w being enforced by the constraint. Given that 
^ /g^° dxx^^2w — vfix^ = 2^ it is straightforward to obtain 



^0 



N' 



N' 



— / dx(j)'^{x) H / dx(j)Q{x)x'^ H ( / dx(j)Q{x)Y 

L 6 ./ 2 ./ 2 ' 



[W W Wl 


= N^w 


[4 + 4 + 2J 







(13) 



as it should, the ground state energy being simply the sum of the zero point 
energies of two (A^^) free oscillators. 

To next order, we assume that the background remains even (J dxx(j)Q = 
0) and obtain from flTTl) 

2A 2A f 

TT^^nfx) = 2fi-2w UJ2 - (w^ H )x'^, 002= I x'^cpoix) (14) 

w w J 

Therefore, to this order, the background remains of the Wigner typo 



7r</)o(x) = V2w — w'^x'^, \x\ < Xo 
It follows that 



2 _2 2 , 2A 
w w 



(15) 



W W Wl 


= N' 


— H h — 


4 4 2J 





w + 



A 



2w^ 



+ ... 



As an example of an equal time correlator, one straightforwardly obtains: 

1 



< TrX'f >= N^ / dxx^6o = N^ — = A^^ 
^ ' ^^ 2w 



_1 X_ 

2^ 2^^'" 



(16) 



These results are in agreement with perturbation theory. 



^fi is suitably adjusted. 



Planar expectation values of correlators involving the X2 coordinate can 
also be calculated. This requires expressing X2 in terms of the B oscillators 
using (jl]) and ([5]), and then reducing the correlator to a correlator involving 
the Xi matrix only, with the use of the B oscillator commutation relations. 
Because the B, B^ sector is quadratic in the oscillators, this can lead to highly 
nontrivial non-perturbative results. For instance. 



<TrXi> = — ^cosh(20,,)-sinh(20y) = i5^ 



1 



2 



^fd.fdy^Mt^^^(l-\^..) (17) 
2 J J ' ^w^ + 2\{x-yY 2w\ -^ ' ^ ^ 



w" 



This equation shows how, to order A, the Xi <-* X2 symmetry of planar 
correlators is satisfied, despite the asymmetric treatment of the two coordi- 
nates in the approach followed in this communication. If we require, as we 
must, that this Xi <->■ X2 symmetry is exact, then fITBl) and flT7|) establish a 
highly non-trivial property of the planar background. 

Of great physical interest is to obtain the states which correspond to the 
spectrum of the system. This can be done in principle using the framework 
developed in [21], [2l] and [25], but is beyond the scope of this communica- 
tion. 

5 Strong coupling solution 

We consider now (ITT]) and fITOl) as A — > 00 : 



TT^cj^lix) = 2/i - 2V2XJ dy\x - ylMv) (18) 



En = N' 



71 f v2A 

dx0o(x) + ^— I dx I dy\x - y\(f)o{x)(f)oiy) (19) 



IQ J ^"^ ' 2 
We find it useful to introduce 



fix) = V2\ I dy\x - y\My), ^Vo(a^) = 2(/i - f{x)) (20) 

which satisfies 



/(^) = ^jdy\^- y|V2(/i-/(y)). (21) 

As it was the case in perturbation theory, we assume that (polx) remains 
an even, single cut function defined in the interval [— a;o,a;o]. To show that 
this is a consistent ansatz, we note that then: 

fix) = V2X{\x\ My)dy + 2 / My)ydy ■ (22) 

Hence f{x) is also even, establishing the consistency of the ansatz. 
Using 

9l\x-y\ = 26{x-y), 
equation (^Th becomecl 

dlfix) = ^v/;i^7M (23) 

n 

This can be integrated in the usual way, to yield: 

kdjr+^-^if^-fix)f^=e (24) 

The "energy" constant is fixed by the condition dxf{0) = 0. Denoting 
/(O) = /o one obtains 



df 4A4 



dx V37r 
The normalization condition 



(/.-io)^-(^-/(x))^ (25) 



fixes 



dx(j)o{x) = 2 / dxcpoix) = 2 / df — ^ — 

■xo Jo Jfo d^ 



(/W-/o)5 = (— )A5, 



t>Q satisfies a very similar equation. 



and hence (1251) takes the form: 



ax 



We will not need to invert (126|) and obtain /(x) explicitly, as all results 
presented in this communication will be expressed in terms of known definite 
integrals. 

Of particular interest is the large A^ ground state energy. From (iT9ll and 
fl20|) this can be written as 



^0 = iV2 



TT- / . ,3, . 1 



— \ rfa;0o(2;) + - / dxf{x)(^Q{x) 



N^ 



L2"T2 Z^^'^ol^)! (27) 



One needs to know /x, or /o, independently. From ( !22l) . one obtains 

/o = V2AX0 - (/U - /o) At = V2Aa;o. 
From (l26l) one obtains 



2Xxo = (/i - /o) 



(i^; 



l-(l-z) 



2(/i - /o 



tdt 



VT^^' 



Also 



TT 
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These integrals are tabulated [27], and are finite. Therefore 



En =N^ 



[?(f)'/ 



tdt 



At 



N' 



^uir-m^'' 



Similarly to the weak coupling case, we consider the correlator 



< TrX^ >=N^ I dxx' 



2N 



= '^xn^df 



/o 



dx 



(28) 



By a sequence of integrations by parts, we obtain 



< TrXl > 



N^ 



2A 



2N^fylx-l 



V2xJfo £ 

1 



N^ 



TC 



2^v^'8 



df 



tdt 



(?)'A-i 



vT^t^ 
3v^/ .2 



vr 



(29) 



rft 



K))°-|(^(^) 



vT^t^ 

3- 



6 Interacting massless matrices 

The results of the previous section could have also been obtained by taking 
the limit w — *> of the planar hamiltonian ([7]). It is therefore relevant to 
discuss the relevance of these results for the system 



H 



1 



Tr{Pl 



+ ^TriR 



gl.,,TT{[X,,X,][X,,X2]), 



(30) 



which results the from the dimensional reduction of massless Higgs or non 
abelian vector potentials. The hamiltonian (130|1 has a single dimensionful 
parameter, A = Qym^ i ^^d therefore all observable quantities should depend 
only on well defined powers of A. However, as is well known, perturbation 
theory is plagued with infrared divergences. In this context, w can thought 
of as a standard "mass" regulator. 

The results of the previous section are therefore remarkable, as they are fi- 
nite and free of any infrared divergences, and depend only on the appropriate 
power of A which is expected from dimension considerations. 

It has already been pointed out the expression ([7]) has a smooth ly — > 
limit. The only place where this limit is potentially ill defined is in the trans- 
formation (jl]), where where w has to be kept finite, if small, to ensure that 
the tranformation is canonical. However, once the Bogoliubov tranformation 
is implemented and the Hamiltonian ([T]) is recast in the form ([2]), its w — *> 
limit should provide a correct description of the system (I30p . 

We are confident that indeed the results of the previous section, which are 
finite and free of any infrared divergences, are the planar energy and corre- 
lator of fl30|) . This provides another explicit confirmation of the expectation 
that strong coupling dynamics, appropriately resummed through the large N 
limit, is free of infrared divergences. 
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7 Summary and discussion 

In this communication, we obtained the large A^ planar background of two 
hermitean matrices, in an harmonic potential, interacting through a Yang- 
Mills potential, in terms of the density of eigenvalues of one of the matrices. 
This background is shown to satisfy a self-consistent gyM dependent integral 
equation. This integral equation has both a weak and a strong coupling 
expansion. The weak coupling expansion of the backround is described and 
shown to be in agreement with perturbation theory. For strong coupling, we 
show that the background satisfies a non-linear differential equation, with a 
solution which is described. The planar ground state energy and examples of 
correlators are obtained and shown to be finite. We argue that this is the full 
planar solution of the hamiltonian of two "massless" matrices (i.e., without 
the harmonic potential, or in the zero curvature limit) with a Yang-Mills 
interaction. 

It may be tempting to associate the background obtained in this com- 
munication with a non-supersymetric qym deformation of the "droplet" de- 
scription of 1/2 BPS states |8], [9], [23]. It should be remembered, though, 
that the matrix description of the geometry is a phase space description. 
The hermitean matrix associated with the droplet corresponds to the holo- 
morphic restriction of a complex matrix [21], unlike the choice made in this 
communicatioEcl 

On the other hand, as pointed out at the end of Section 4, by exploiting 
the Xi <->■ X2 symmetry of the system, we can write from (fT7|l and as A — >■ 00, 

< TrX^ >=< TrXl >= - J] ^ • (31) 

With a physical interpretation of the eigenvalues as coordinates of a system 
of DO's, the feature of fl29|) that its "size" becomes small but finite for large 
A can be interpreted as resulting from the "back-reaction" exhibited on the 
right hand side of ( l3Ti) . 

Clearly, it would be very interesting to examine this question further. 
However, at the very least, the background identified in this article is a qym 
deformation of the Wigner distribution associated with harmonic potentials, 
with a well defined strong coupling limit. 



^As an example of how the kinetic energy operator of the haniihonian is changed in 
this case, see l28l. 
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